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In water management, it is essential to know the water availability 
within the hydrological planning and management units. To improve the 
knowledge of hydrological cycle variables in the Zirahuén Lake Basin we 
used the Soil-Water-Balance (SWB) spatiotemporal simulation model 
supported by the hydrodynamic characteristics of the soil and land 
cover, together with daily climatological information. Our results showed 
how land use influences runoff and recharge of aquifers reflecting in 
water availability estimates, which showed an important setback with an 
emphasis during the 2009-2015 period. When comparing our results 
with the official records reported in 2016, we found that the latter 
figures underestimated the groundwater index (GI) by 14.64% and the 
relative water availability index (RWAI) by 150% relative to the results 
of our study. We concluded that, to improve the reliability and 
management of water in the basin, it is recommend verifying our results 
by implementing strategies to monitor the quantity, quality, and 
temporality of water. 


Keywords: Water Management, Runoff, Groundwater, Water balance, 
Water use. 


Resumen 


En la gestión del agua es fundamental conocer la disponibilidad hídrica 
al interior de las unidades de planeación y manejo hidrológico. En la 
cuenca del lago de Zirahuén, para mejorar el conocimiento de las 
variables del ciclo hidrológico, se empleó el modelo de simulación 
espaciotemporal en cuencas Soil-Water-Balance (SWB), apoyado en las 
características hidrodinámicas del suelo y cobertura, en conjunto con 
información climatológica diaria. Los resultados mostraron la influencia 
del uso del suelo en los escurrimientos y la recarga de acuíferos, que se 
reflejaron en la disponibilidad hídrica, la cual mostró retrocesos 
importantes con énfasis en el periodo 2009-2015. Al analizar los 
resultados con respecto a lo reportado en 2016, se encontró que el 
Indice de Aguas Subterráneas (IDAS) ha sido subestimado en 14.64%, 
en tanto que la Disponibilidad Relativa (DR) se ha subestimado por el 
150% con respecto a los valores encontrados en este estudio. 
Finalmente, para mejorar la confiabilidad y la gestión del agua, se 
recomienda implementar estrategias de monitoreo de la cantidad, 
calidad y temporalidad del agua para verificar los resultados obtenidos. 
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Introduction 


Water and its availability are the base for development of societies, 
although its distribution lacks a relationship with population growth and 
expansion, and with economic development (Shiklomanov, 1998). In 
Mexico, renewable water availability is in an inverse proportion with 
socioeconomic level ( Conagua, 2016a)), which generates high-demand 
zones having no trustworthy information regarding the functioning of 
the water cycle. 


The water balance (WB) is the result of the evaluation of the factors 
involved in the water cycle, for which methods have been developed for 
different observation scales (Sokolov 8 Chapman, 1974; Xu € Singh, 
1998). The soil-water balance (SWB) model ( Dripps 8: Bradbury, 2007; 
Westenbroek, Kelson, Dripps, Hunt, €: Bradbury, 2010) estimates the 
components of the water balance at a daily level (Thornthwaite 8 
Mather, 1957) taking advantage of geographic information systems. It 
calculates the water balance at the basin scale using five required 
variables: 1) Precipitation and temperature, 2) soil use and vegetation, 
3) hydrological soil group, 4) flow direction, and 5) water retention 
capability of soil. The runoff values are calculated by means of the 
runoff curve numbers of the National Resources Conservation Service 
(NRCS, 1986). The method has been successfully used for regional 
groundwater recharge evaluation in the Lake Michigan and Wisconsin 
basins (Tillman, Gangopadhyay, €: Pruitt, 2016). 


The joint analysis of the water balance components and the involved 
water demand allow to know water availability (Touhami et a/., 2014), 
which in Mexico is estimated by the Mexican official norm NOM-011- 
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Conagua-2015 (Conagua, 2015) and is only existing from 2014 to the 
present with lack of data for several time periods. In order to improve 
the management of water and due to the lack of information through 
time, the simulation models allow improving the knowledge of the water 
balance (and its availability) in basins. Because of that, the objectives of 
our work were to apply the SWB model for estimating water availability 
in the Zirahuén Lake Basin at different times, and to analyze the 
availability indexes relative to the reported official records for 2016. 


Materials and methods 


The study was made using the geospatial information of the Zirahuén 
Lake Basin (19%21'09" and 19%29'50" N, and 101%29'11" and 
101%46'21" W) in the state of Michoacán de Ocampo (Figure 1) within 
the Tarasca Neovolcanic Physiographic Province (INEGI, 2001). 
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Figure 1. Location of the Zirahuén Lake Basin. 


Information was taken from the 1:50,000 scale topographic map 
(INEGI, 2014) to generate a digital elevation model (DEM) at a 100 m 
spatial resolution. In the process of information conditioning, the land 
cover and land use classes (Table 1) were converted to ASCII raster 
format data matrices. A database from The National Commission of 
Water (Conagua) with daily information of temperature and precipitation 
was used converting values to English units of measurement. 


Table 1. Synthesis of the used information. 


Type Format Version Scale 
Topographic Vector Series III 1:50,000 
Soil Vector Series II 1:250,000 
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Basin limits Vector 2nd edition 1:50,000 
Soil use and soil cover Vector Series 1, 11, 11, | 1:250,000 
IV, V 
Climatological data Table 2016 N/A 


The evaluation of the spatiotemporal variations in the water balance was 
made for the years 1972, 1980, 1992, 2000, 2003, 2009, and 2015 
based on the publication date of each land use and land cover series and 
the availability of daily climatological data. The evapotranspiration was 
estimated using data of maximum and minimum temperatures 
(Hargreaves 8 Samani, 1985) taking into  consideration the 
meteorological information in the daily records. 


To calculate the water availability (surface and groundwater), the Water 
Rights Public Record (REPDA) database was used (Conagua, 216c) that 
only contains data of discharge volumes from the year 2000 onwards, 
because of which the calculations were made for a similar period. For 
the natural discharge of groundwater, the same value (0.3 Mm?) of the 
“Nueva Italia” polygon was used because the basin's area is within that 
polygon. Regarding the return volume of surface water, we assumed a 
conservative criterion considering that the lowest value found (2015) for 
the return volume (1.45 Mm?) was constant. 


The groundwater index (GI) was calculated based on the results in the 
REPDA database and for the relative water availability index (RWAI) we 
used the data reported by Conagua (2016b) in the basin substituting the 
natural runoff value with the data found in this work. 


In general, the GI and RWAI values were calculated as established in 
the NOM-011-Conagua-2015 for the year 2015 and compared with the 
official results reported for the beginning of 2016 (Conagua, 2016b). 


The land use and land cover classes were coded based on the reference 
values table of the SWB model (Westenbroek et al., 2010) to be 
converted to raster and to generate maps of land use change dynamics 
during the analyzed period (Figure 2). 
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Figure 2. Land cover and land use coded for the SWB. 


Four soil groups were identified in the Series II 1:250,000 scale soil 
units and profile points (INEGI, 2004): Andosol, Leptosol, Luvisol, and 
Phaeozem. Three textural classes were detected by crossing the 
information from the soil units and soil profiles, one for each 
hydrological group (NRCS, 1986): Andosol with a sandy-loam texture in 
hydrological group A with a code value of 1, Leptosol with a clayey- 
sandy-loam texture in group C and code value 3, and Luvisol and 
Phaeozem with a clayey-loam texture in group D code 4. 


The values of the components of soil texture were used to assign values 
based on the soil type tables (Thornthwaite € Mather, 1957). A 1.6 
value was assigned to Andosols, 2.7 to Leptosols, and 3 to Luvisols and 
Phaeozem. 
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Results 


The recharge depth per year (Figure 3), showed the highest values for 
the years 1992 and 2000 and the lowest recharge values for the years 
2009 and 2015. The north-central part covered by forests had low 
recharge vales due to the steep slope present in the higher elevations. 
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Figure 3. Spatial and temporal distribution of recharge. 
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The precipitation values in 1972 were similar to those in 2015, but the 
recharge values were in the order of three-fold lower than in 1972 
(Table 2), displaying the effects of land use change. 


Table 2. Components of the water balance in recharge depth per year 


NC-SA 4.0 
c-sa/4.0/) 


(mm/y). 
Year | Precipitation Recharge Interception | Runoff Potential 
evapotranspiration 
(PET) 
1972 999.00 178.82 26.92 484.12 1 238.25 
1980 1 219.60 156.46 27.18 596.90 1771.65 
1992 1922.10 415.80 32.51 935.23 1 560.07 
2000 1 406.00 206.76 31.24 681.74 1 695.70 
2003 1 071.80 107.70 44.45 519.68 1555.50 
2009 958.89 53.09 43.18 465.58 1613.15 
2015 999.50 45.21 41.40 487.43 1 636.27 


The water availability in the basin (Table 3) is the result of the detailed 
analysis of the REPDA database. A 284% increase was observed in the 
demand of groundwater from 2009 to 2015 motivated by the volume of 
water required for agriculture, which in that period went from 0.1033 
Mm? to 0.3975 Mm”. Regarding the demand of surface water, the public 
supply experienced a 40% increase between 2009 and 2015. These 
results show there is an intensification of the pressure on water within 
the basin. 


Table 3. Water availability in millions of cubic meters per year. 


283 
Tecnología y ciencias del agua, 10(3), 275-288. DO!: 10.24850/j-tyca-2019-03-11 


2019, Instituto Mexicano de Tecnología del Agua 
7 Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y 


CienciaszAgua (https://creativecommons.org/licenses/by-nc-sa/4.0/) 
Year R Di DHi R D2 DH2 

2000 56.52 0.10 56.12 186.36 0.75 187.06 

2003 29.44 0.19 28.95 142.06 2.19 141.32 

2009 14.51 0.18 14.03 127.27 2.44 126.29 

2015 12.36 0.75 11.31 133.25 3.18 131.52 


Abbreviations: R, recharge; D1, demand on groundwater extraction; 
DH1, groundwater availability; R, runoff; D2, demand on surface 
extraction; DH2, surface water availability. 


Discussion 


We explained the observed differences in the recharge depth per year 
by the soil use changes (Figure 2) occurred in the forest, grasslands, 
and cropland areas, which affected the recharge rates. It has been 
found that land use changes from forest to agriculture have an effect on 
soil and water conservation, as well as causing a decrease in the 
infiltration rates (Sun et al., 2018). In degraded basins, crop rotation 
applied as a corrective measure reduces soil erosion by up to 69% 
(Abdelwahab, Bingner, Milillo, 8 Gentile, 2014); otherwise, the 
conversion from forest cover to croplands (agriculture and foraging) 
causes a depletion of the organic carbon in mineral soils and of carbon 
in rootlets (Mendoza-Vega, Karltun, € Olsson, 2003). Land use change 
in the scenario of climate change decreases the annual runoff discharge 
by 19% and the maximum summer flow by up to 45% (Pijl, Brauer, 
Sofia, Teuling, € Tarolli, 2018). 


The 'Nueva Italia” aquifer has a Gl value of 0.82 and is classified within 
the availability zone 4, with enough volume to satisfy the demand of 
water occurred during 2016 (Conagua, 2016b). The calculations based 
on the simulation (2015) resulted in Gl values of 0.94, 14.16% more 
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than the official estimates for the total aquifer. Analyzing the granted 
water volumes for the 2000-2015 time period (Table 3), there was a 
linear trend with a negative slope that indicates a decrease along time in 
the amount of available water from the aquifer. 


Regarding surface water, the results from the simulation displayed a 
decreasing trend with a minimum in 2009 despite the apparent recovery 
occurred during the year 2015; the values were less than those for 2000 
and 2003. The official figures show a RWAI value of 0.92 within the 
availability zone 1; substituting the runoff volumes resulting from the 
SWB model (Table 3) the obtained RWAI was 2.32, which is over 150% 
more than the official calculations, thus placing it within availability zone 
2, in which there is enough available volume to satisfy the present 
needs and a margin for future demand. 


Our calculated runoff values (133.246 Mm? in 2015) are far larger than 
the estimations made by Conagua (52.16 Mm”) in the determination of 
water availability, the capacity of the basin having been 
underestimated; therefore, the water availability is 78.755 Mm?, These 
differences are mainly due to the method used to estimate the runoff: 
Conagua uses the rational method with annual results, while the SWB 
model uses the runoff curve numbers in a geospatial solution taking into 
account the particular topography and land cover class, obtaining daily 
water balance results. The differences in the time lapse of the 
calculations allow a higher sensitivity of the results. Yet, in both 
processes, the information about the amount of water is obtained 
without considering the quality of water. 


Although we used a more detailed modeling, the scope of our work has 
not been verified by gauging because, until today, there is a void of 
information in that regard. Likewise, the information about the quality 
and temporal distribution of water is insufficient. The SWB model has 
proven its usefulness in the simulation of the water cycle with relevant 
results in the U.S.A. (Tillman et a/., 2016; Zhang, Felzer, 8 Troy, 2016), 
having a variation of +33 mm per year in the estimation of the 
recharge, and a high confidence level has been found in the results of 
the model at the regional scale when using data at a monthly and yearly 
resolution (Han et al., 2018). The studies about the application of the 
SWB model account for its usefulness as a tool to support the processes 
of water management within basins. 
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Conclusions 


The use of the SWB model allowed estimating the water balance with a 
good approximation of the groundwater and surface water flows. As a 
result of taking into consideration the biophysical characteristics of the 
basin and the precipitation at a daily temporal resolution, the analysis 
regarding water availability allowed adjusting the GI values 14.92% 
above the official figures and the RWAI values at over 150% higher than 
the official reports. Although the literature reports the usefulness of the 
SWB model in the estimation of recharge, our results show the urgent 
need to implement gauging instrumentation in the Zirahuén Lake Basin, 
and also of strategies for the monitoring of the quantity, quality, and 
temporality of water in order to improve the process of its management. 
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